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Abstract 

A method for obtaining the small current quark mass effect on the dressed quark propagator 
from an effective quark-quark interaction model is developed. Within this approach both explicit 
and dynamical chiral symmetry breaking has been analyzed. A comparison with previous results 
is given. 
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The dynamical chiral symmetry breaking in QCD is a typical nonperturbative phe- 
nomenon, where the standard perturbative scheme is not valid. While lattice QCD is the 
most straightforward and solid approach among nonperturbative studies of QCD, analytical 
studies with effective theories are useful alternatives and often shed more light on nonper- 
turbative phenomena than lattice QCD. The calculation of the last twenty years show that 
the global color symmetry model(GCM) provides a successful description of various non- 
perturbative aspect of strong interaction physics[l-4]. We naturally expect that it might be 
a useful model in the study of the small current quark mass effects on the dressed quark 
propagator which manifests the behavior of the explicit and dynamical chiral symmetry 
breaking. 

The previous study of dressed quark propagator for the light quarks u and d has been 
done often in the chiral limit, where the current quark mass is set to zero. However, in a more 
realistic situation, i.e., quarks have small but finite masses representing explicit symmetry 
breaking which play an important role in the study of QCD phase structure. It is interesting 
to address small current quark mass effects on the dressed quark propagator further. Up 
to this let us start from the Euclidean action of the GCM at small but finite current quark 
mass m: 



SGCM[q,q;m] 

— j d'^x {q{x)[y ■ dx + m]q{x)} + j d'^xd'^y 

where j'^{x) = q{x)^^^q[x) denotes the color octet vector current and g'^D'^Kx — y) is the 
dressed model gluon propagator in GCM. Without loss of generality, we will employ the 
Landau gauge D^^{q) — 5ab {^^v — D{q^) for the model gluon propagator. 

Introducing an auxiliary bilocal field B^{x,y) and applying the standard bosonization 
procedure the partition function of GCM [1,2] 

Z[m] = J pgpge"^'^^^ (2) 

can be rewritten in terms of the bilocal fields B^{x, y) 

Z[m] = J VB^e-^^ff^'''''^^ (3) 

with the effective bosonic action 

Seff[B^; m] = -Trlng~^[B'^; m] + /[S^] (4) 



and the quark operator 

g-^[B'; m] = [7 • a, + m]5{x - y) + E(x, y), (5) 

where the exphcit expressions of I[B^] and = A^B^{x,y) can be found in Ref.[5]. 

The matrices = D"^ ® ® F'^ is determined by Fierz transformation in Dirac, color and 
flavor space of the current current interaction in Eq.(l). 

In the mean-field approximation, the fields B^{x,y) are substituted simply by their vac- 
uum value Bq{x, y), which is defined as 



= 0. (6) 



5B 

Employing the stationary condition Eq.(6), and reversing the Fierz transformation, we have 
quark self-energy Eo(x, y) at the mean-field approximation 

4 

Eo[m\{x,y) = -glD^^{x,y)jf,go[m]{x,y)j^. (7) 

It should be noted that both B^^x^y) and Go^{x^y) are dependent on the small current 
quark mass m. If m is switched off, ^o["^] goes into the dressed vacuum quark propagator 
in the chiral limit G = Qo[m = Q], which has the decomposition 

G-\p) =ij-p + So(p) = ^7 ■ pMp') + Bip") (8) 

with 

So(p) = / d'xe^^-[K' Bl{x)] = y ^d^AP - QhMlh^ (9) 

where the self energy functions A[p^) and B{jp) are determined by the rainbow Dyson- 
Schwinger equation(DSE) in the chiral hmit: 



W)-lV-lj-^dD{p-^) 



^p- {p- q) q- {p- q) 
p-q + 2 

[P - q) 



(27r)4^^ IP 5J ^2^2(^2) +52(^2)- 

Here we should stress that the B{p'^) in Eq.(lO) has two qualitatively distinct solutions. 
The "Nambu-Goldstone" solution, for which B{p^) ^ 0, describes a phase, in which: 1) chiral 
symmetry is dynamically broken, because one has a nonzero quark mass function; and 2) the 
dressed quarks are confined, because the propagator described by these functions does not 
have a Lehmann representation. The other solution, the "Wigner" one, B{p^) = 0, describes 



a phase, in which chiral symmetry is not broken. In "Wigner" phase, the Dyson-Schwinger 
equation (10) reduces to: 



3 J (27r; 



p ■ q + 2 



p- {p-q) q- (p-q) 



(11) 



{p - qf 

where A^p^) denotes the dressed quark vector self energy function in the "Wigner" phase. 
Therefore, in chiral limit, the dressed quark propagator in the "Wigner" phase can be written 
G^'^Kq) = 

Let us now study the small current quark mass dependence of the dressed quark prop- 
agator. To this, for the light quarks u and d one can expands in powers of m as 
follows 



+ 



which leads to the formal expansion 



m + 0(m^)^G-^ + mr + 0(m^), (12) 



m=0 



gJm] =G-mGTG + 



with r 



r(i/i,2/2) = 



Sm 



(13) 



(14) 



m=0 



In coordinate space the dressed vertex r{x, y) is given as the derivative of the inverse quark 
propagator ^(^^[m] with respect to the m. 

Taking the derivative in Eq.(5) and putting it into Eq.(12), we have 

5So[m] (7/1,7/2) 



r(7/i, ^2) = S{yi - y2) + 



5m 



(15) 



m=0 



Substituting Eq.(7) and (13) into Eq.(15), we have the inhomogeneous ladder Bethe- 
Salpeter equation(BSE) of scalar vertex F, which reads 

r (1/1,2/2) =5(2/1-2/2) 

4 f 

-^9sDnuiyi-y2) J duidu2'yi^G{yi,ui)r{ui,U2)G{u2,y2hu- (16) 
Fourier transform of Eq.(16) leads then to the momentum space form of F 

F(p, 0) = 1 - ^ / -0^yMp - qhMqMq, ^)G{q)l.. (17) 

As was shown above, both the rainbow DSE(IO) and the ladder BSE(17) can be consistently 
derived from the action of the GCM(l) at the mean field level. 



Prom Lorentz structure, the most general form for the F which fulfills Eq.(17), reads: 

r(p,0) = z7-pC(p2) + D(p2)_ (18) 

Substituting Eq.(18) into Eq.(17), we have the coupled integral equation for C(p^) and 



^glD{p-q) 



p - q + 



2p-{p-q) q-{p-q) 
{p-qY 



X 



/ [2A{q^)B{q^)D{q') + q'A\q')C{q') - C{q^)B\q^ 
I [q'A^q^) + B^{q^)f 

d^q 



Dip') -l-4j^y,D(p-q)x 



[2q^A{q^)B{q^)C{q^) - q^A\q^)D{q^) + D{q^)B^{q^)] 



(19) 



1 [gM2(g2) + 52(^2)]2 J- 

In "Wigner" phase, B{p^) = 0. Substituting B{p^) — into Eq.(19), we have the 
decoupled integral equation for C'{p^) and D'{p^) in "Wigner" phase: 



1 + 4 



d^q 



glD{p-q)- 



p- q + 

■ D'{q' 



2p-{p-q) q-{p-q) 
(p - qf 



C'{q') 



(20) 



(27^)4^^^^/^ ^^^2^/2(^2-)- 

Prom Eq.(20), it is easy to find C'ip^) = in the "Wigner" phase. 

Based on the above coupled integral equation(19) and Eq.(lO), by means of numerical 
studies, we can get the nonperturbativc scalar vertex F, which is useful for the studies of 
the small current quark mass effects on the dressed quark propagator. It should be noted 
that the above approach for getting the nonperturbativc scalar vertex F has been proven 
to be very useful for the studies of nonperturbativc vector, axial vector vertex[6-8] and the 
chemical potential dependence of the dressed quark propagator [9]. So far, at the mean field 
level, we have completed the derivation of the first order dependence of upon m in 

"Nambu-Goldstone" and "Wigner" phase separately; 



= i'y-p ^(p^) + 5(p^) + m fi7 • p C(p^) + £'(p^)l + C(m^), 



(21) 



^ i^-p A'{p'^) + mD'{p'^) + 0{m'^). 



(22) 



= i'j ■ p A(p^) + B(p^) is the contribution of dynamical symmetry breaking 
and is independent of the current quark mass m. ^q^*^^ (E) = m [^7 ■ p CijP') + -D(p^)] is 
the contribution of the exphcit symmetry breaking and vanishes if the current quark mass 
m = 0. It is now apparent that the above approach for getting the nonperturbative scalar 
vertex has the advantage that the effects of explicit and dynamical chiral breaking can be 
analyzed separately. 

As is shown in Eq.(21), for m 7^ the dressed-quark self energies in general acquire not 
only the scalar part D{p^) but also the vector part C(p^) driven by the current quark mass 
m. One often miss the contribution of the vector part C(]?^) in the previous study of dressed 
quark propagator for light quarks. 

In order to have a qualitative understanding of the contribution of C(p^), -D(p^), and 
D'(p^), a particularly simple and useful model of the dressed gluon two-point function[10] 
is used as follows: 



g'iD^^ip - g) = Any 



S^'\p-q). 



(23) 



(p-q)^ 

Here the scale parameter 1] is a measure of the strength of the infrared slavery effect. It 
should be noted that the model gluon propagator(23) is an infrared-dominant model that 
does not represent well the behavior of ggDf^i,{p) away from ~ 0. 

Then substituting Eq.(23) into Eq.(lO) and (11), we have the "Nambu-Goldstone" solu- 
tion; 
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S(p') = (?7^-V)2, Aip')^2 far p' < 



B{p' 



0, A{p^ 



1+1 + -^^ 



for 



4 ' 

2 ^ 
P >j, 



and the "Wigner" solution; 



B'ip^) = 0, A'ip^) = \ 



1+1 + ^ 



(24) 



(25) 



With the model of the dressed gluon propagator specified in Eq.(23) and Eqs. (24,25), 
from Eq.(19) and (20) we have the following simple form of the C{p'^), D{p^) 



Dip') = 



p^ + if +p^{l + 2if/p^)^ 

p2 _ ^2 -\-p^{\ + 27]^ Ip"^)^ 



P'>\, 



(26) 



lr^2 + 



2 772 - 4p2 ' 



C{p') = 



2{ri^ -Ap^y. 
rf — 4p2 



P2<^ 



4 ' 



and C'(f), D'(p') 



D'ip^) = l±t±fil±lilf}i, CV)-0. (27) 



2 



At large momentum region, i.e., — > oo, from Eq.(26), we have 

Dip") 1, Cip") = 0. (28) 

This asymptotic behavior is quite different from that given by Ref.[ff]. This is because 
the model gluon propagator specified by Eq.(23) is an infrared-dominant model that does 
not represent well the real large momentum behavior of glDni,{p) in QCD. This question 
will be further discussed by means of a more sophisticated ansatze for the gluon propagator 
elsewhere. 

With these two "phase" characterized by qualitatively different momentum-dependent 
quark propagator(21,22), the GCM can be used to calculate the vacuum condensates and 
susceptibilities at the mean field level[12]. In particular we obtain the chiral quark conden- 
sate with the small current quark mass m: 

(0| : qq : |0)^ = -tmc {^^f^^ N - Q^'''^ [m] ] (29) 

3 nV^ r 2(ry2 - 4s) \2{rf - Asf/"^ + m{Tf + 8s)] 

ass 



4tt^ Jo [ 16s[(7y2 - 4s) - m{rj^ - 4s)V2]2 + [2(^2 _ 4^)3/2 + ^(^2 + g^^^ 

2m[s -rf + s(l + 2rf / sf/"^] | 
~ [s - 772 + s(l + 2r;7s)V2]2 + 2m?[s + r]^ + s(l + 2r;Vs)V2] j ' 

Substituting m = into Eq.(29), we have the chiral quark condensate in chiral limit: 

(0| : qq : |0)^=o = -tvnc {^f ""^[m = 0]} = (30) 

The calculated ratio (0| : qq : |0)^/(0| : qq : 10)^=0 can be seen from Fiq.l. In Fig.l, we 
see that (0| : qq : |0),„/(0| : qq : |0),„=o increases with increasing m within the certain range 
of the small current quark mass (it should be noted that our numerical result is only valid for 
low values of m). In order to show the contribution of the subtraction term troc {Oo^\ni\^ , 
we plot the ratio troc {^o^''["^]} / (0| • QQ '■ |6)m=o versus m in Fig.2. It is easy to find that 
the contribution of subtraction term can not be neglected in calculating the chiral quark 
condensate with finite current quark mass m. 
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Fig.l. The ratio (0| : qq : |0)m/(0| : qq : |0)m=o as a function of m for r]"^ = 1.14 GeV^. 
I ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 
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Fig. 2. The ratio tr^c {^o^'*["^]} / (0| : qq '■ |0)m=o as a function of m for rj'^ = 1.14 GeV'^. 

To summarize: in order to discuss exphcit and dynamical chiral symmetry breaking, 
we first provide a general recipe to calculate the small current quark mass effect on the 
dressed quark propagator at the mean field level in the framework of GCM. This employs a 
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consistent treatment of dressed quark propagator G and the dressed scalar vertex F, which 
are both determined from the effective quark-quark interaction by the rainbow DSE for G 
and the inhomogeneous ladder BSE for F. This approach has the advantage that we can 
analyze the effects of explicit and dynamical chiral symmetry breaking separately. Then 
we use a simple, confining model(23) to calculate the chiral quark condenste in the case 
of the small current quark mass. It is found that the contribution of the subtraction term 
tvDC {po^\^]\ can not be dropped. 
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